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Abstract 

A search for double /? decay of dysprosium was realized for the first time with the help 
of an ultra low-background HP Ge 7 detector. After 2512 h of data taking with a 322 
g sample of dysprosium oxide limits on double beta processes in ^^^Dy and ^^^Dy have 
been established on the level of T112 > 10^^ — 10^® yr. Possible resonant double electron 
captures in ^^^Dy and ^^^Dy were restricted on a similar level. As a by-product of the 
experiment we have measured the radioactive contamination of the Dy203 sample and 
set limits on the a decay of dysprosium isotopes to the excited levels of daughter nuclei 
as ri/2 > 10^5 - 10^^ yr. 



PACS: 23.40.-s; 23.60. +e 
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1 INTRODUCTION 



The double beta (2/3) decay experiments are considered to-date as the best way to determine 
an absolute scale of the neutrino mass and to establish the neutrino mass hierarchy, to clarify 
the nature of the neutrino (Majorana or Dirac particle), to look for existence of right-handed 
admixtures in the weak interaction and of hypothetical Nambu-Goldstone bosons (Majorons), 
and to test some other effects beyond the Standard Model [1]. The developments in the new 
experimental techniques during the last two decades lead to an impressive improvement of 
sensitivity to the neutrinoless (Oz/) mode of 2/3~ decay up to 10^^ — 10^^ yr [2]. Allowed in the 
Standard Model the two neutrino (2i/) double beta decay was detected for 10 nuclides with the 
half-lives in the range of 10^^ - 10^^ yr [21 [3]. 

^Corresponding author. E-mail address: rita.bernabei@roma2.infn.it (R. Bernabei). 
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The sensitivity of the experiments to search for the double electron capture (2e), the electron 
capture with emission of positron {e(3~^), and the double positron (2/3+) decay is substantially 
lower: the best counting experiments give only limits on the level of 10^^ — 10^^ yr [21 El |6p. 
There is a strong motivation to develop experimental technique to search for these processes: 
study of neutrinoless 2e and e(3~^ decays could clarify the contribution of the right-handed 
admixtures in weak interactions [8]. Dysprosium contains two potentially 2/3 active isotopes: 
i56j3y -vyii;]^ one of the largest releases Q213 = (2012 ± 6) keV (therefore both 2e and eP~^ 
channels of decay are possible), and ^^^Dy {Q213 = (284.6 ± 2.5) keV, only double electron 
capture is energetically allowed) [9]. The decay schemes of the triplets ^^^Dy — ^^^Tb — ^^^Gd 
and ^^®Dy — ^^^Tb — ^^^Gd are presented in Fig. 1 and 2, respectively. 
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Figure 1: Decay scheme of ^^^Dy [lO}. Energies of excited levels and emitted 7 quanta are in keV 
(relative intensities of 7 quanta, roughed to percent, are given in parentheses). Excited states 
related with possible population and subsequent de-excitation of 0+ levels at Eexc = 1851.2 
and 1715.2 keV, and levels at E^xc = 1827.8 and 1771.1 keV are not shown. Deexcitation 
branches with low probabilities (< 1%) also are not shown. 



■^An indication for 2/3+ decay processes in ^■^"Ba and ^^^Ba was obtained in geochemical measurements [7]; 
however, this result has to be confirmed in a direct counting experiment. 
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Figure 2: Decay scheme of ^^^Dy |TT]. Energy of excited levels and 7 quanta are in keV. 



It should be mentioned the possibility of a resonant enhancement of the neutrinoless double 
electron capture in ^^^Dy and ^^^Dy due to energy degeneracy. The resonant double electron 
capture was discussed in Refs. [121 [131 EH ES] ) where an enhancement of the decay rate by some 
orders of magnitude was predicted for the case of coincidence between the released energy and 
the energy of an excited state. According to [15], high Z atoms are strongly favored to search 
for resonant 2e decay. Dysprosium has one of the highest Z among nuclides for which resonant 
processes could occur. 

Resonant captures are possible on a few excited levels of ^^^Gd and one level of ^^®Gd. The 
properties of the excited levels are listed in Table [H Because transitions with difference in spin 

Table 1: Characteristics of possible resonant neutrinoless double electron capture in ^^^Dy and 
^^®Dy. Difference — -Eexc (where Q2/3 is the double beta decay energy, ii^exc is the energy of 
an excited level) is denoted as 5. 



Transition 


Decay channel 


Energy (keV) and parity 
of excited level 


S (keV) 


156Dy^l56Gd 


IK 


1914.8, 2+ 


-3 ±6 




KLi 


1946.4, 1- 


7±6 




KLi 


1952.4, 0- 


1 ±6 




2Li 


1988.5, 0+ 


7±6 




2L3 


2003.8, 2+ 


-6 ±6 


158l3y_^158Gd 


2Li 


261.5, 4+ 


6.3 ±2.5 



more than 2 are strongly suppressed, we consider in this study only the levels of ^^^Gd with 
spin < 2. However, we left in the list the level of ^^^Gd with the spin 4"*" to which the resonant 
capture is possible. 

Unfortunately, the isotopic abundances of both potentially double beta active dysprosium 
isotopes are rather low: concentrations of ^^^Dy and ^^^Dy in the natural dysprosium are 
0.056(3)% and 0.095(3)%, respectively [IS]. 

To our knowledge there were no attempts yet to search for double (3 decays of ^^^Dy and 
^^^Dy. The aim of the present work was the search for 2/3 processes in the dysprosium isotopes 
with the help of ultra-low background high purity (HP) Ge 7 spectrometry. As a by-product 
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of the experiment we have estimated the radioactive contamination of the dysprosium oxide 
sample and set limits on the a decay of the dysprosium isotopes to excited levels of the daughter 
nuclei. 

2 MEASUREMENTS, RESULTS AND DISCUSSION 

2.1 Contamination of Dy203 measured by mass-spectrometry 

A sample of high pure dysprosium oxide (Dy203 of 99.98% purity grade) was provided by 
the Institute for Single Crystals (Kharkiv, Ukraine). The contamination of the material was 
measured with the help of High Resolution Inductively Coupled Plasma-Mass Spectrometric 
analysis (HR-ICP-MS, Thermo Fisher Scientific ELEMENT2). Potassium was measured in 
High Resolution mode. The calibration of the apparatus was performed by adding in one of 
the sample solutions of a known amount of a standard solution containing this element. All the 
other elements were determined in Low Resolution High Sensitivity mode. A semiquantitative 

Table 2: Contamination of Dy203 analyzed by ICP-MS analysis. The calculated activities of 
the most important radioactive isotopes presented in the radioactive elements are given. 



Element 


Concentrations 


Activity of 




of element 


radioactive isotopes 




(ppm) 


(mBq/kg) 


K 


< 0.5 


< 15 (40K) 


Ca 


50 




Co 


< 0.05 




Cu 


< 5 




Sr 


< 0.02 




Y 


6 




La 


0.25 


0.21 (i38La) 


Sm 


17 


2200 (i^^Sm) 


Gd 


5 


0.008 (i52Gd) 


Tb 


< 5 




Ho 


30 




Er 


60 




Yb 


3 




Lu 


0.5 


26 (i^^Lu) 


Pb 


0.16 




Th 


0.036 


147 (232Th) 


U 


0.004 


2.3 (235U), 49 (238U) 



analysis was performed, that is a single standard solution containing some elements at a known 
concentration level (10 ppb of Li, Ce, Tl) was used for calibration and quantification. In this 
case, the uncertainty can be estimated as about 30% of the given concentration value. The 
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results of the HR-ICP-MS analysis are presented in Table [H We have also calculated the 
activity of the radionuclides present in the measured elements. 

2.2 Low background measurements 

A Dy203 sample with the mass of 322 g in a polyethylene bag was placed directly on the end 
cap of the ultra-low background HP Ge detector (GeBer, 244 cm^) installed deep underground 
(~3600 meters of water equivalent) at the Gran Sasso National Laboratories of the INFN 
(Italy). The detector was located inside a passive shield made of low radioactivity lead (~ 20 
cm), copper 10 cm) and borated polyethylene (~ 10 cm). To remove radon, the set-up 
was continuously flushed by highly pure nitrogen. The energy resolution of the spectrometer 
is 2 keV at 1332 keV 7 line of ^^Co. The Dy203 sample was measured over 2512 h, and the 
background of the detector was collected during 6110 h. 

The energy spectrum accumulated with the Dy203 sample by the HP Ge detector is shown 
in Fig. 3 together with the background data. The spectra are normalized on the times of 
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Figure 3: (Color online) Energy spectra accumulated with dysprosium sample over 2512 h 
(Dy203) and without sample over 6110 h (BG) by ultra-low background HP Ge 7 spectrometer. 
(Inset) Low energy part of the spectra. Energy of 7 lines are in keV. 

measurements. The peaks in the spectra can be ascribed to 7 quanta of the radionuclides ^°K, 
^°Co, and daughters of U/Th. We have also observed 7 lines of ^''^Lu (201.8 and 306.8 keV) 
in the data accumulated with the dysprosium sample. To estimate radioactive contamination 
of the Dy203 sample, the detection efficiencies of the background components were simulated 
with the help of the EGS4 package [17j. The initial kinematics of the particles emitted in the 
decay of nuclei was given by an event generator DECAYO [18]. Radioactive contamination of 
the dysprosium oxide is presented in Table El The data agree with the results of the ICP-MS 



5 



Table 3: Radioactive contamination of Dy203 sample measured by HP Ge 7 spectrometry. 



Chain 


Nuclide 


Activity (mBq/kg) 






< 10 




6OC0 


< 1 




137Cs 


< 2 




i38La 


< 3 




i^^Lu 


9±2 


232T11 




179 ± 40 




^^^Tli 


157 ± 15 


235 -jj 


235 u 


7±3 






< 80 




227AC 


< 190 


238 u 


234mpg^ 


< 115 






5±3 




210pb 


< 9400 



analysis taking into account accuracy of the measurements. 

The observed contamination of the sample by U/Th is quite predictable taking into account 
that the principal ores for production of dysprosium are monazite and bastnaesite. Monazite 
contains considerable (typically a few %) amount of thorium and uranium. The presence of 
radioactive ^^^Lu can also be expected due to the very similar chemical properties of lutetium 
and dysprosium, which lead to certain difficulties in chemical separation of the elements. 

2.3 Search for 2e and e/S^ decay of ^^^Dy 

There are no clear peculiarities in the energy spectra accumulated with the dysprosium sample, 
which can be interpreted as double beta decay of ^^^Dy or ^^^Dy. Therefore only lower half-life 
limits can be set according to the formula: 



lim Ti/2 = N ■ 1] ■ t ■ ln2/ lim S, 

where is the number of potentially 2/3 unstable nuclei, rj is the detection efficiency, t is the 
measuring time, and limS* is the number of events of the effect searched for which can be 
excluded at a given confidence level (C.L.). 

To estimate values of lim S for 2f3 processes in ^^^Dy and ^^^Dy, the energy spectrum ac- 
cumulated with the Dy203 sample was fitted in different energy regions where 7 peaks related 
with the decay processes are expected. The efficiencies of the HP Ge detector for the double (3 
processes searched for were simulated with the help of the EGS4 code [T7] and the DECAYO 
event generator [18]. 

One positron can be emitted in the ef3~^ decay of ^^^Dy with the maximal energy of (990 ±6) 
keV. The annihilation of the positron will give two 511 keV 7's leading to extra counting rate 
in the annihilation peak. The energy spectra accumulated with and without the sample in the 
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Figure 4: (Color online) Part of energy spectra accumulated with Dy203 sample over 2512 h 
by ultra-low background HP Ge 7 spectrometer (Dy203). The energy spectrum accumulated 
without sample over 6110 h normalized on the time of measurements with the Dy203 sample 
(BG). Fits of 511 annihilation 7 peaks are shown by solid lines. 



energy interval (400 — 700) keV are presented in Fig. 4. The measured area of the annihilation 
peak is (862 ±44) counts. The area of the peak in the background spectrum is (149 ±23) counts 
(normalized on the time of measurement with the Dy203 sample). The excess of events in the 
data accumulated with the dysprosium sample can be explained by radioactive contamination 
of the material. Indeed decays of ^^^Ra and ^^^Th daughters give in total (757 ± 72) counts 
in the annihilation peak. The difference between the observed and the expected number of 
counts (—44 ± 87) gives no evidence for the effect of the decay of ^^^Dy. In accordance 
with the Feldman- Cousins procedure [I9], we obtain 102 counts which can be excluded at 90% 
C.L. (here and below all the half-life limits and values of limS* are given at 90% C.L.). Taking 
into account the number of ^^^Dy nuclei in the sample {N = 2.91 x 10^°) and the detection 
efficiency (3.4%), we have obtained the following limit on the half- life of ^^^Dy relatively to the 
sum of 2i/ and Oz/ modes of e(3~^ decay: 

y(2.+o.)e/3+^i56j3y^g_g_ ^ g.s.) > 1.9 X lO^*^ yr. 

The e(3~^ decay of ^^^Dy is also allowed to the first excited level of ^^^Gd with the energy of 
89.0 keV. To estimate an upper limit on the transition, the energy spectrum was fitted in the 
energy interval (83 — 109) keV by the model constructed from a Gaussian function at the energy 
of 89 keV with the energy resolution FWHM = 1.6 keV (the 7 peak searched for), a linear 
function describing the background, and a second Gaussian to take into account the 7 peaks 
with the energies 92.4 and 92.8 keV (^'^^Th from ^ssu chain). A fit by the chisquare method 
{x^/n.d.f. = 18.9/19 = 0.99, where n.d.f. is number of degrees of freedom) gives the area of 
the peak searched for: S = (86 ± 88) counts, which gives no evidence for the effect. We took 
230 counts which can be excluded at 90% C.L. Taking into account the simulated efficiency 
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to 89 keV 7 quanta (0.05%) we have obtained the hmit on 2z/£/3"'" decay of ^^^Dy to the first 
2+ 89.0 keV excited level of ^^^Gd: Tllf^{g.s. 89.0 keV) > 1.3 x 10^^ yr. However, a 
much stronger restriction can be set by considering the annihilation peak. Indeed, the detection 
efficiency for the 511 keV 7 for the 2i'e[3^ (and QveP^) decay to the 89.0 keV excited level is 
3.4% which leads to the same limits as for the e^^ decay to the ground state: 

y(2.+o.)./3+^i56j3y^ g.s. ^ 89.0 keV) > 1.9 x lO^^ yr. 

In a case of the 2v2K capture in ^^^Dy a cascade of X rays and Auger electrons with the 
individual energies up to 50 keV is expected. The most intensive X ray lines are 42.3 keV 
(26.6%), 43.0 keV (47.5%), 48.6 keV (4.8%), 48.7 keV (9.3%) and 50.0 keV (3.1%) |20]. To 
derive a limit on the two neutrino double electron capture from K shells, the energy spectrum 
was fitted in the energy interval (35 — 60) keV by the model consisting of five Gaussians (to 
describe the expected X ray peaks) plus two Gaussians (46.5 keV 7 line of ^^°Pb and 53.2 
keV 7 line of ^^^U), and a polynomial function of the first degree (background; see Fig. 5). 
The fit gives the total area of the 2v2K effect as (310 ± 102) counts. It is worth noting 
that the background in this energy region is quite complex to be safely estimated; therefore, 
conservatively we give here only a limit on the number of events: lim 5* = 477 counts. Taking 
into account the detection efficiency of the effect (0.5%), one can calculate the following half-life 
limit: 

Tyi^{^^^T)Y, g.s. ^ g.s.) > 6.1 X 10^^ yr. 

To estimate a limit on the 2v2£ decay of ^^^Dy to the excited levels of ^^^Gd, the energy 
spectrum accumulated with the dysprosium sample was fitted in energy intervals where peaks 
from the de-excitation 7 rays are expected. For instance, we have obtained the following limit 
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Figure 5: (Color online) Low energy part of energy spectrum accumulated with Dy203 sample 
over 2512 h. An excluded effect of 2v2K process in ^^^Dy (^^^Dy) with the half-life T1/2 = 
6.1 X 10^^ yr (T1/2 = 1-0 x 10^^ yr), and 2v2e transition of ^^^Dy to the first excited level of 
^^^Gd (89.0 keV, T1/2 = 1.8 x 10^^ yr) are shown by sohd lines. 

on the 2v2e decay of ^^^Dy to the first excited 89.0 keV level of ^^^Gd: 
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T^/f ^(i56j3y^ g.s. ^ 89.0 keV) > 1.8 x 10^^ yr. 

A peak expected for the 2z/2£ transition of ^^^Dy to the excited 2+ 89.0 keV level of ^^^Gd 
is shown in Fig. 5. 

The limits obtained for the 2i/2e decay of ^^^Dy to the excited levels of ^^^Gd are presented 
in Table m where the energies of 7 quanta used in the analysis, detection efficiencies and values 
of lim^* are given too. 

In case of 0z/2e decay to the ground state of ^^^Gd, we suppose here that only one bremsstrah- 
lung 7 quantum is emitted to carry out the transition energy (in addition to X rays and Auger 
electrons from de-excitation of atomic shell). The energy of the 7 quantum is expected to be 
equal to = Q213 — E^i — Eb2, where £"61 and Eb2 are the binding energies of the first and of 
the second captured electron of the atomic shell. The binding energies on the K and Li, L2, 
shells in gadolinium atom are equal to E^ = 50.2, E^^ = 8.4, E^^ = 7.9 and E^^ = 7.2 keV 
[20] . respectively. Therefore, the expected energies of the 7 quanta for the 0^26 capture in 
^s^Dy to the ground state of ^^^Gd are: i)E^ = (1912 ± 6) keV for 0i^2K; n) E^ = (1954 ± 7) 
keV for OuKL; ni) E^ = (1996 ± 7) keV for 0i/2L. 

There are no clear peaks in the energy intervals (1906 — 1918), (1947 — 1959), and (1989 — 
2004) keV expected in the g.s. — )■ g.s. 0^26 decay of ^^^Dy (see Fig. 6). We have estimated 
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Figure 6: (Color online) Part of the energy spectra where peaks from 0i/2e processes {2K, 
KL and 2L) in ^^^Dy to the ground state of ^^^Gd are expected. The excluded peaks for the 
processes are shown. 



values of lim by a fit of the data in the energy intervals by the simple model: a Gaussian 
function (to describe the peaks searched for) plus polynomial function (continuous background). 
Taking into account the calculated efficiencies to detect 7 quanta with energies ~ 1.9 — 2.0 MeV 
(~ 0.7%), we set the following limits on the processes searched for: 



rp0u2K (156-^^^ g-s. ^ g.s.) > 1.7 X 10^^ yr, 
jOuKL(i56-Q^^ g.s. ^ g.s.) > 1.7 X 10^6 yr, 

g.s.) > 2.2 X 10^^ yr. 



y0.2L(156j3y^ g.s. 
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The high Q213 energy of ^^^Dy allows also population of several excited levels of ^^^Gd with 
energies in the range of (89 — 2004) keV. To estimate limits on these processes, the energy 
spectrum accumulated with the dysprosium sample was analyzed in different energy intervals. 
A summary of the obtained results is given in Table HI 

2.4 2£ decay of ^^^Dy 

An estimation of the half-life of ^^^Dy relatively to the two neutrino double K electron capture 
was obtained by using the limS* and the detection efficiency for the 2v2K decay of ^^^Dy. 
Taking into account the number of ^^^Dy nuclei in the sample (A^ = 4.94 x 10^°) we have set 
the following limit on the 2u2K decay of ^^^Dy: 

rp2u2K ^ibsj^^^ g-s. g.s.) > 1.0 X 10^5 yr. 

The simulated response of the detector to the 0i/2£ decay of ^^^Dy to the ground state of 
^^^Gd differs from the two neutrino process. We have set the following limit by analyzing the 
data in the vicinity of the 184.2 keV 7 peak expected for the 0z/2A' decay: 

^vf^l^^^Dy, g-s. ^ g-s.) > 4.2 X 10^6 yr. 

Limits on the 2v and Oz/ modes of double electron capture to the first excited level 79.5 keV 
of ^^®Gd were obtained by analysis of the data in the energy region around the energy 79.5 keV 
(see Table i]). 

2.5 Search for resonant double electron capture in ^^^Dy and ^^^Dy 

The double electron capture of ^^^Dy to the excited levels of ^^^Gd with energies of 1914.8, 
1946.4, 1952.4, 1988.5 and 2003.8 keV may occur with higher probability due to a possible 
resonant enhancement. The limits on the resonant double electron capture in ^^^Dy from K 
and L shells were obtained by analyzing the experimental spectrum in the energy intervals 
where the most intensive 7 peaks from de-excitation of these levels are expected. For instance, 
we set the following limit on the resonant 2 A' capture in ^^^Dy to the excited level 2"*" of ^^^Gd 
with the energy 1914.8 keV (the dysprosium spectrum was fitted in the vicinity of the 7 line 
1826.0 keV): 

y^(2^^+o;.)2£j.i56j3y^ g.s. ^ 1914.8 keV) > 1.1 x 10^^ yr. 

To search for a resonant double electron capture in ^^^Dy from two Li shells to the excited 
level 261.5 keV of ^^®Gd, the experimental spectrum was fitted in the energy interval 175 — 185 
keV. The fit gives an area of a 181.9 keV peak (—48 ± 49) which corresponds to limS* = 40 
counts. The detection efficiency for 7 quanta with the energy 181.9 keV irradiated in the 2e 
decay is 1.3%, therefore we can set the following limit on the process: 

^(2.+o.)2e^i58j3y^ g.s. ^ 261.5 keV) > 3.2 x lO^^ yr. 
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Table 4: Half-life limits on 2^ processes in ^^^Dy and ^^^Dy. 



Process 


Decay 


Level of 


Ej 


Detection 


lim^* 


Experimental 


of decay 


mode 


daughter 


(keV) 


efficiency 




limit (yr) 






nucleus 














(keV) 


























2z/ + 0z/ 


g.s. 


511 


3.4% 


102 


> 1.9 X 10^6 




2i/ + 0i/ 


2+ 89.0 


511 


3.4% 


102 


> 1.9 X 10^*^ 


2K 


2y 


g.s. 


42.3 - 50.0 


0.5% 


477 


> 6.1 X 10^^ 


2e 


2u 


2+ 89.0 


89.0 


0.07% 


230 


> 1.8 X 10^^ 


2e 


2v 


0+ 1049.5 


960.5 


1.1% 


9 


> 7.1 X 10^6 


2e 


2v 


2+ 1129.4 


1040.5 


0.6% 


24 


> 1.4 X 10^^ 


2e 


2v 


2+ 1154.1 


1154.1 


0.5% 


62 


> 4.7 X 10^5 


2e 


2v 


0+ 1168.2 


1079.2 


1.0% 


65 


> 8.9 X 10^5 


2e 


2v 


0+ 1715.2 


472.7 


1.1% 


21 


> 3.0 X 10^6 


2e 


2v 


2+ 1771.1 


1682.2 


0.7% 


39 


> 1.0 X 10^6 


2e 


2v 


2+ 1827.8 


1738.9 


0.2% 


6 


> 1.9 X 10^6 


2e 


2v 


0+ 1851.2 


697.0 


0.1% 


39 


> 1.5 X 10^5 


2K 


Oi/ 


g.s. 


1906 - 1918 


0.7% 


24 


> 1.7 X 10^6 


KL 


Oi/ 


g.s. 


1947 - 1959 


0.7% 


24 


> 1.7 X 10^6 


2L 


Oz/ 


g.s. 


1989 - 2001 


0.7% 


18 


> 2.2 X 10^6 


2e 


Oz/ 


2+ 89.0 


89.0 


0.06% 


230 


> 1.5 X 10^^ 


2e 


Oi/ 


0+ 1049.5 


960.5 


1.0% 


9 


> 6.4 X 10^6 


2e 


Oi/ 


2+ 1129.4 


1040.5 


0.6% 


24 


> 1.4 X 10^6 


2e 


Oi/ 


2+ 1154.1 


1065.2 


0.5% 


70 


> 4.1 X 10^5 


2e 


Oz/ 


0+ 1168.2 


1079.2 


0.9% 


65 


> 8.0 X 10^^ 


2e 


Oz/ 


0+ 1715.2 


472.7 


1.0% 


21 


> 2.8 X 10^6 


2e 


Oi/ 


2+ 1771.1 


1682.2 


0.6% 


39 


> 8.9 X 10^^ 


2e 


Oi/ 


2+ 1827.8 


1738.9 


0.2% 


6 


> 1.9 X 10^6 


2e 


Oi/ 


0+ 1851.2 


697.0 


0.1% 


39 


> 1.5 X 10^^ 


Resonant 2K 


2z/ + 0i/ 


2+ 1914.8 


1826.0 


0.5% 


27 


> 1.1 X 10^6 


Resonant KLi 


2z/ + 0z/ 


1- 1946.4 


1857.4 


0.4% 


24 


> 9.6 X 10^^ 


Resonant KLi 


2i/ + 0z/ 


0^ 1952.4 


709.9 


1.2% 


27 


> 2.6 X 10^6 


Resonant 2Li 


2i/ + Oi/ 


0+ 1988.5 


1899.5 


0.7% 


21 


> 1.9 X 10^6 


Resonant 2L3 


2i/ 


2+ 2003.8 


1715.1 


0.02% 


41 


> 2.8 X 10^^ 


Resonant 2L3 


Oi/ 


2+ 2003.8 


761.3 


0.03% 


58 


> 3.0 X 10^^ 
















2K 


2v 


g.s. 


42.3 - 50.0 


0.5% 


477 


> 1.0 X 10^5 


2K 


Oi/ 


g.s. 


184.2 


1.7% 


40 


> 4.2 X 10^6 


2e 


2z/ 


2+ 79.5 


79.5 


0.04% 


113 


> 3.5 X 10^^ 


2e 


Oz/ 


2+ 79.5 


79.5 


0.03% 


113 


> 2.6 X 10^^ 


Resonant 2Li 


2i/ + 0i/ 


4+ 261.5 


181.9 


1.3% 


40 


> 3.2 X 10^6 
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All the half-life limits on 2/3 decay processes in dysprosium obtained in the present experi- 
ment are summarized in Table HI 

It should be stressed that the difference between the decay energy and the energy of an 
excited level is one of the most crucial characteristics for the resonant double electron capture. 
However, in most of the cases the accuracy of atomic mass data is still not enough to point 
out isotopes with minimal differences Q2/3 — E^xc- Therefore searches for resonant transitions 
in different isotopes are strongly motivated. 



2.6 Search for a decay of dysprosium isotopes 

Five of the seven natural dysprosium isotopes are potentially unstable for a decay. Furthermore 
for all the isotopes transitions to excited levels of daughter nuclei are allowed. Characteristics 
of the possible alpha decays are presented in Table [51 The search for a decay of the dysprosium 
isotopes to the excited levels of the daughter gadolinium isotopes was realized by analysing the 
low background energy spectrum measured with the Dy203 sample. 

We do not observe any peculiarities in the experimental data which could be interpreted as 
7 peaks from a decays of the dysprosium isotopes to the lowest excited levels of the daughter 
nuclei, and set only limits on these decays. To estimate a limit on the a decay of ^^^Dy to the 
excited 2"^ 344.3 keV level of ^^^Gd, the energy spectrum of the Dy203 sample was fitted in 
the energy interval (341 — 350) keV by a simple model: Gaussian function (to describe 7 peak 
with the energy of 344.3 keV) and polynomial function of the second degree (to approximate 
background, see Fig. 7). A fit by the chisquare method (x^/n.d.f. = 5.3/5 = 1.1) gives 
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Figure 7: (Color online) A part of the energy spectrum of HP Ge spectrometer accumulated 
with Dy203 sample over 2512 h. An expected 7 peak from a decay of ^^^Dy to the first excited 
level of ^^^Gd corresponding to the half-life T1/2 = 4.3 x 10^^ yr (excluded at 90% C.L.) is shown 
by the solid line. 



S = (—48 ± 42) counts for the area of the peak searched for, which provides no indication for 
the effect. According to the Feldman- Cousins procedure [19] we took 31 counts as lim S. Taking 
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into account the detection efficiency (2.1%) and the yield of 7 quanta from the transitions from 
the excited level 2+ 344.3 keV (0.962), we have obtained the following limit on the a decay of 
^^^Dy to the first excited level of ^^^Gd: 

Tf/2(^^''Dy, g.s. ^ 344.3 keV) > 3.8 x 10^^ yr. 

Limits on T1/2 relatively to a decay to the lowest excited levels of the daughter nuclei for 
other dysprosium isotopes were obtained in the same way. The results are presented in Table [5] 
where theoretical estimations calculated by using the approach proposed in [21] are given too. 

Table 5: Characteristics of possible alpha decays of dysprosium isotopes to the lowest excited 
levels of the daughter nuclei. 



All 1 

Alpha decay 


T— > • J_ 1 1 1 

bxcited level 


Detection 


1 * 

lim b 


bxperimental 


rni j_ • 1 

1 heoretical 


Qa (keV) [9] 


01 daughter 


efficiency 




1 * * J 

limit on 


estimation of 




nuclei (keV), 
yield of 7 

per a decay 


of the 7 
quanta 




J- 1/2 lyrj 


m 




2+, 344.3 


2.1% 


31 


> 3.8 X 10^6 


2.3 X 10^^ 


1759(6) 


0.962 |2I] 










0.056(3) 












158l3y^l54Gd 


2+, 123.1 


0.9% 


31 


> 1.3 X 10^6 


1.4 X lO^s 


879(3) 


0.455 |22] 










0.095(3) 












160j3y^l56Gd 


2+, 89.0 


0.4% 


230 


> 8.5 X 10^^ 


5.4 X 10^26 


439.2(13) 


0.203 [in] 










2.329(18) 












16lDy^l57Gd 


5/2-, 54.5 


0.2% 


83 


> 3.5 X 10^6 


1.3 X 10^4^ 


344.6(12) 


0.075 [23] 










18.889(42) 












162j3y^l58Gd 


2+, 79.5 


0.3% 


113 


> 1.0 X 10^^ 


1.7 X 10^^14 


85.0(12) 


0.142 [11] 










25.475(36) 













The theoretical half-lives were found very big. Indeed, the predicted half-life is T1/2 = 2.3 x 10^ 
yr even for the decay with the biggest energy release (^^^Dy — )■ ^^^Gd, 2"*", 344.3 keV). 

Theoretical estimates give no hope for realistic observation of alpha decays of all naturally 
abundant dysprosium isotopes to the excited level of daughter nuclei. At the same time the 
detection of the ^^^Dy decay to the ground state of ^^^Gd (supposing one will apply crystal 
scintillators from dysprosium enriched in ^^^Dy) looks possible. Indeed we have estimated the 
half-life as 5 x 10^^ yr by using the approach [2l|. The important advantage of scintillation 
methods is an almost 100% detection efficiency to alpha particles and pulse-shape discrimination 
ability [251 126]. Development of enrichment methods for dysprosium isotopes and radiopure 
crystal scintillators containing Dy are requested to observe a activity of natural dysprosium. 
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3 CONCLUSIONS 



The first experiment to search for 2/3 processes in ^^^Dy and ^^^Dy was carried out in the 
underground Gran Sasso National Laboratories of INFN by using ultra-low background HP Ge 
7 spectrometry. After 2512 h of data taking with a 322 g sample of highly pure dysprosium 
oxide (Dy203) limits on double beta processes in ^^^Dy and ^^^Dy have been established on the 
level of 10^^ - 10^6 yr. 

The search for the resonant 0^26 capture in ^^^Dy and ^^^Dy to excited levels of daughter 
nuclei was of particular interest. Taking into account the strong dependence of the process on 
the difference between the double beta decay energy and the energy of an excited level, a more 
accurate determination of the difference in atomic masses of ^^^Dy - ^^^Gd and ^^'^Dy - ^^^Gd 
isotopes is strongly required. A precise study of the 1914.8, 1946.4, 1952.4, 1988.5 and 2003.8 
keV levels of ^^^Gd characteristics (spin, parity, decay scheme) is important too. 

Theoretical calculations of double beta decay of dysprosium are scarce; in fact, the only 
known estimations for ^^^Dy (g.s. to g.s. transitions) are T^j2^ = 7.0 x 10^'' yr (for the 

effective neutrino mass of 1 eV) [27J, T^f^ = (1.1 - 2.8) x 10^6 yr, T^^^ = (1.8 - 4.5) x 10^3 
yr [2B] • The present study is the first experimental attempt to search for double beta processes 
in ^^^Dy and ^^^Dy, which both are rather interesting isotopes for 2/3 experiments taking into 
account a possible resonant double electron capture on excited levels of daughter nuclei. 

The sensitivity of the experiment can be advanced to the level of T1/2 ~ 10^'^ — 10^^ yr 
by using enriched ^^^Dy and ^^^Dy isotopes deeply purified from radioactive elements, and 
increasing the exposure and detection efficiency by application of multi-crystal HP Ge detectors. 

The search for alpha decay of dysprosium isotopes to the excited levels of daughter nuclei 
was realized as a by-product of the experiment on the level of sensitivity limTi/2 ~ 10^^~^^ 
yr. However, the obtained limits are very far from the theoretical predictions even for the 
a transition of ^^^Dy with the highest energy of decay (T1/2 ~ 10^^ yr). At the same time 
the detection of a decay of ^^^Dy to the ground state of ^^^Gd (theoretical estimations of the 
half-life are on the level of ~ 10^^ yr) looks reasonably realistic by using crystal scintillators 
containing dysprosium enriched in ^^^Dy. Such crystal scintillators could be also used to search 
for double beta processes in ^^^Dy. 

We have found that the Dy203 sample contains on the level of several mBq/kg ^^^Lu, 
^3^U and ^^^Ra; the activities of ^^^Ra and ^^^Th are ^ 0.18 and f» 0.16 Bq/kg, respectively. 
The results of the 7 spectrometry are consistent with the data of the ICP-MS analysis for 
radionuclides. The contamination of the dysprosium sample can be explained by the production 
of lanthanides from monazite, an ore with high concentration of uranium and thorium. Presence 
of radioactive ^''^Lu can be due to the similar chemical properties of lutetium and dysprosium, 
which provides some difficulties in chemical separation of the elements. 
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